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ABSTRACT 

We analyzed the host stars of the present sample of confirmed planets detected by Kepler and Kepler 
Objects of Interest (KOI) to compute new photometric rotation periods and to study the behavior 
of their angular momentum. Lomb-Scargle periodograms and wavelet maps were computed for 3, 807 
stars. For 540 of these stars, we were able to detect rotational modulation of the light curves at 
a significance level of greater than 99%. For 63 of these 540 stars, no rotation measurements were 
previously available in the literature. According to the published masses and evolutionary tracks of the 
stars in this sample, the sample is composed of M- to F-type stars (with masses of 0.48-1.53 Mq) with 
rotation periods that span a range of 2 to 89 days. These periods exhibit an excellent agreement with 
previously reported (for the stars for which such values are available), and the observed rotational 
period distribution strongly agrees with theoretical predictions. Furthermore, for the 540 sources 
considered here, the stellar angular momentum provides an important test of Kraft’s relation based 
on the photometric rotation periods. Finally, this study directly contributes in a direct approach to 
our understanding of how angular momentum is distributed between the host star and its (detected) 
planetary system; the role of angular momentum exchange in such systems is an unavoidable piece of 
the stellar rotation puzzle. 

Keywords: — stars:— general — Hertzsprung-Russell and C-M diagrams — stars: kinematics and 
dynamics — stars: rotation 


1. INTRODUCTION 

The study of stellar rotation is important for our un¬ 
derstanding of many kinematic properties of stars, in¬ 
cluding the origin and evolution of the rotation itself; 
more generally, the origin and evolution of stellar an¬ 
gular momentum; and the role of rotation in stellar 
magnetism and mixing of chemical elements. Stars are 
formed from rotating molecular cloud cores and pre¬ 
serve only a very small fracti on of the initial angular 
momentum of these cores (e.g., lPallal[2?)?)^ iLamm et al.l 
1200511 . Indeed, the initial angular momentum of the 
initial molecular cloud is approximately eight orders of 
magnitude greater than that of the stars that eventu¬ 
ally form from such a cloud core (e.g., iBodenheimerl 
I989H. With the disco very of extra-solar planets (e.g., 
Mayor fc Qnelod 1 199511 . new questions concerning rota¬ 
tion have begun to arise, including questions regard¬ 
ing the large-scale effects of th e presence of planetary 
companions ( 3n ste llar rotation (jPontl 120091 : lAlves et all 
l20irt iLanzal [201011 and on t he heating of stellar coro- 
nae and chromospheres (e.g.. iPoppenhaeger et al.ll2010l 

iCanto Martins et al.ll201lt ILanzal 201211 . Admittedly, the 

study of the rotation of planet-hosting stars is ham¬ 
pered by one major factor: for the great majority of 
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stars with planets only the projected rotational veloc¬ 
ity V sini is available, and all such bodies identified, 
to date are essent ially slow-to-very-moderate rotators 
(jAlves et al.l 1201(1( 1. Despite this constraint, Alves et 
al. have noted that stars without any detected planets 
demonstrate a clear angular momentum deficit compared 

with stars with planets. __ 

The recent space-missions C oRoT ((Baglin et al.ll2007tl 
and Kepler (iKoch et al.ll201(ill are providing a new win¬ 
dow into the study of the behavior of stellar rotation, 
including the rotation of stars with planets, by virtue 
of their high-precision photometry, which yields excel¬ 
lent quality light curves (LCs). In addition to their 
high precision, these missions provide nearly uninter¬ 
rupted photometric measurements of unprecedented du¬ 
ration and cadence, allowing for refined analysis of the 
features in volved in the behavior o f these temporal se- 
ries, as in Ide Freitas et ah (l2013b[l. In recent works , 
McQuillan et al.l (l2013allbl . l2014^Nielsen et al.l (I20I3I1 
Walkowicz fc Basril (1201311 . and iReinhold et al.l (l2013f l 
computed the rotation periods (Prot) for a large fraction 
of the Kepler targets, i ncluding the Kepler Obje cts of In- 
terest (hereafter KOI). iMcQuillan et al.l (|2013a]lbl . I20I4H 
analyzed the first 11 quarters of Kepler data by comput¬ 
ing Fourier periodograms and autocorrelation functions 
(AC F), and validated their res ults through visual inspec¬ 
tion. IMcQuillan et al.l (|2013ari computed the rotation pe¬ 
riods of 1, 570 M dwarfs and identified a possible change 
in the slo pe of the period-mass rela tion at approximately 
0.55 Mq. IMcQuillan et al.l ()2013br i measured the periods 
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of 737 main-sequence KOI stars to study the relation be¬ 
tween the rotational and orbital periods; they found that 
hotter stars rotate more rapidly. iMcQuillan et al.l (|2014ll 
derived periods for 34 030 main-sequence stars, exclud¬ 
ing eclipsing binaries and KOI stars, and found that 
the majority of these stars are typically younger than 
4.5 Gyr and exhibit higher pea k-to-peak am p htude s 
for smaller pe ri ods. More over. iNielsen et al.l (l2013f) . 
iReinhold et al.l (|2013ll . and iWalkowicz fc Baso (1201^ 
computed pe riods automat i cally u sing Lomb-Scargle pe- 
riodograms. INielsen et al.l (1201311 report e d rot ation pe¬ 
riods for 12,151 stars. IReinhold et al.l (j2013l l derived 
periods for 24,124 targets using an automatic method 
with partial visual inspection. These authors identi¬ 
fied a trend of magnetic braking when comparing rota- 
tion periods w ith the gyrochronological relations from 
iBarnesI (j2007[ l . IWalkowicz fc Basril (120131 1 reported rota¬ 
tion periods and ages for approximately 950 KOI stars, 
and through comparison with spectroscopic data, they 
estimated the corresponding stellar orientation angles. 
These authors also used gyrochronological ages based on 
BarnesI (1200711 . in combination with the calibrations of 
Mamaiek fc Hillenbrandl ()2008l l to estimate stellar ages. 

In this work we present rotation periods for 131 planet¬ 
hosting Kepler stars. We computed the rotation periods 
using a robust procedure, t hat combine s the standard 
Lomb-Scargle periodogram (iLombI 119761: iScargle 1198^ 
with the wavelet technique (jGrossmann fc Morletl 198411 . 
We also present rotation periods for a comparison sam¬ 
ple of 409 KOI stars (released December 13, 2013). By 
combining these data with theoretical predictions, we are 
able to present an overview of the rotational behavior of 
these stars on the Hertzsprung-Russell (HR) diagram. 
This work is organized as follows: in Section^ the work¬ 
ing stellar sample, the Kepler data and the procedure 
used for the computation of periods are presented. Sec- 
tion|3]presents the results. Our conclusions are presented 
in Sectional 


2. KEPLER DATA AND ANALYSIS 

From May 2009 to May 2013 the NASA Kepler mis¬ 
sion collected data in an steady field of view for 191,449 
stars in 17 runs (known as quarters) which were com¬ 
posed of long cade nce (6.02 s observati ons stacked ev¬ 
ery 29.4 minutes, iJenkins et al.l l2010b) and short ca¬ 
dence (bins of 59 secorids) o bservations (I Van Cleve et al.l 
120101 : [Tompson et al.ll201^ . We selected all public light 
curves LGs available in the Kepler Public Archival for 
planet-hosting stars and KOIs between quarters 1 and 
17, which were collected over the noted range of four 
years (from May 13, 2009, to May 08, 2013). From this 
parent sample, 408 Kepler Confirmed Planetary Host 
Star^ (hereafter KCP) were available in the Kepler Pub¬ 
lic Archive. Details regarding the p ublic data have been 
discussed in many publica. t ions (e.g..lBorucki et al.ll2009L 
2010 : iBatalha et al.l[20l^ iKoch et al.l 120101: iBasri et^ll 
201 Ih . In particular, the Kepler database provides 


both Simple Aperture Photometry (SAP) data, which 
were processed using a standard treatment, and Pre- 
Search Data Conditioning (PDC) data, which were an¬ 
alyzed using a more refined treatment based on the 


^ http://archive.stsci.edu/kepler/ 

® http://kepler.nasa.gov/Mission/discoveries/ 
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Figure 1. Example of a LC before (upper panel) and after (bot¬ 
tom panel) the removal of transits and thermal noise. 

PDC routine from the Kepler pipeline (| Jenkins et alJ 
l2010ajl . The PDC routine primarily removes thermal 
and kinematics effects cau sed by the spacecraft (see 
iVan Cleve fc Caldwellll200^ for further details). 

In this work, the parent samples of the KCP host stars 
and KOIs were analyzed using the PDC LCs. The KOIs 
were used as a comparison sample, because of the large 
number of stars with this classification (almost all KCP 
host stars are also KOIs). Most of these LCs are of ex¬ 
cellent quality and required no additional treatment of 
artifacts. However, for the final solutions some flare- 
lik^ signatures (as described by iMaehara et al.l 1201211 
and known planetary transits were removed. Removing 
transits and flare-like signatures allows for better isola¬ 
tion of other sources of stellar LC variability, which may 
include rotational modulation. In total we removed all 
transits from all 408 stars of the KCP parent sample and 
61 flare-like signatures found in 26 (of the 408) stars of 
the KCP parent sample. We also removed outliers by dis¬ 
carding all flux measurements that exceeded 3.5 cr from a 
S^d-order polynomial fit to each quarter’s LC. Further¬ 
more, so me few remaining instrum ental trends (as de¬ 
scribed bv iPetigura fc Marcvll2012ll were removed. These 
included, in particular, thermal noise (which was pre¬ 
dominantly produced by variations in the overall space¬ 
craft temperature that were induced by the monthly 
pointing of the antenna toward Earth for data trans¬ 
mission), which manifested as exponential decays in th e 
LCs (see iPetigura fc Marc^ 120121 : iGarcfa et all 1201 ID . 
and safe-mode signatures (induced by cosmic rays hit¬ 
ting t he detector), whic h also exhibited high-slope de¬ 
cays (|Carcfa et al.ll20lTl l. Both types of artifact exhib¬ 
ited very similar shapes, and they were visually identified 
for exclusion. Figure [T] presents an example of a LC from 
which transits and artifacts were removed. Because both 
the Lomb-Scargle periodograms and wavelet maps were 
calculated using discrete transforms of the time series, 
both types of analysis were capable of treating the gaps 
introduced by removing the data, as described above. 
Finally, individual quarters were combined to produce 
a single long-term LC for each object. The quarter- 
to-quarter transitions typically manifested as flux off- 

® We named flare-like signature a sudden and strong flux bump 
in the LC in short time interval (typically of few days), whose 
physical or instrumental origin is not important to this work. 
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sets (jumps), which were adjusted through a linear fit 
and extrapolation of user-def ined boxes before and af - 
ter each jump, as describ ed bv iDe Medeiros et alJ (|2013[1 
and iBanvai et alJ (|2013r i. 

For data analysis, Lomb-Scargle periodograms were 
computed over the entire s ample within a frequ ency 
range of 0.01 - 0.5 d“^. As in iRdnhold et al.l (1201311 . we 
defined a maximum frequency of 0.5 d“^ (or a minimum 
period of 2.0 d) to avoid pollution of the samples with 
sources of variability other than rotational modulation 
(typically pulsations or oscillations). T he Lomb-Scargle 
metho d was based on the approach of iPress fc Rvbica 
(1198911 and standard statistica l methods fe.g.. iScargld 
119821 : iHorne fc Baliuna.^ 119861 : IPress fc Rvbickil 11989(1 
were used to compute confidence levels. Accordingly, 
in each periodogram, we identified the main peak with 
a confidence level of greater than 99%. We also com¬ 
puted wavelet maps and global wavelet power spectra 
within the same frequency range that was used for the 
Lomb-Sc argle periodograms , following the procedure de¬ 
scribed in iBravo et al.l ^201411 . These authors treated the 
6‘^-order Morlet wavelet as the mother wavelet, because 
of its good time localizatio n and frequency resolution 
(|Grossmann fc Morletl HOS^l . The width of the mother 
wavelet was variable and was defined by Si = 
(|Torrence fc Comooll 19981 1. where Si is the so-called scale 
(the period decomposition in our case), sq is the mini¬ 
mum scale (with a value of 2.0 for our calculations), 5 is 
the scale bin (with a value of 0.1 for our calculations), 
and rii is an iterative index that ranges from 0 to A^ — 1, 
where N is the total number of scales. Wavelet anal¬ 
ysis is a powerful tool for analysis of frequency varia¬ 
tions over time, for a given signal decomposed at var¬ 
ious resolutions. In addition to revealing the stability 
of superposed variability periods over time, the time- 
frequency representation of wavelet maps allows us to 
identify specific signatures that are related to variability 
behaviors. Typical examples of the wavelet maps com¬ 
puted in our analysis are displayed in Fig. [51 together 
with their corresponding LCs. For the statistical analy¬ 
sis of global wavelet periodograms, we followe d the ap¬ 
proach presented in iTorrence fc Comool (119981 1 (see § 4 
and Eq. 20 in the cited publication). We then identified 
the main peaks in the global wavelet periodograms with 
confidence levels of greater than 99%. 

Subsequently, we performed an automatic pre¬ 
selection of KOIs and KCP host stars, which were then 
visually inspected for our final selection. This automatic 
pre-selection was less strict for KCP host stars than for 
KOIs because the KCP host stars comprised our main 
sample, which was sufficiently small for visual inspection 
to be feasible. Below we describe the procedures used 
in our pre-selection for the KCP and KOI subsamples. 
Finally, we estimated the stellar angular momentum of 
each target (the results are presented in § [3]) and also 
estimated the uncertainties of all computed parameters. 
These final steps are described below. The computed pe¬ 
riods, with values that ranged from 2.64 - 36.25 d (KCP 
host stars), and 2.05 - 88.94 d (KOIs), are available in 
Tabled! 


2.1. Pre-seleetion of KCP 


For our main sample, harmonic fits of the LCs were 
computed to estimate their va r iabilit y amplitudes, as de¬ 
scribed in iDe Medeiros et al.l (|2013fl . and the fit resid¬ 
ual was assumed to be the LC noise. Thus, the 
peak-to-peak amplitude over the standard deviation of 
the LC noise was defi ned as the variabi l ity si gnal-to- 
noise S/N ra tio fsee iDe Medeiros et"^ 1201311 . Sim- 
ulations from [Pe Medeiros et al.l ( 201311 (see Fig. 2 of 
IDe Medeiros et al.ll2013ll suggest that a reasonable pre¬ 
selection can be obtained by using a cutoff ot S/N > 1.0, 
when the sample is sufficiently small to be visually in¬ 
spected. 

After this automatic pre-selection, the resulting sub¬ 
sample was visually inspected by apply i ng th e list of 
criteria described by IDe Medeiros et ^ (|2013ll . which 
serve to define the photometric signature of rotational 
modulation caused by dynamic star spots (namely semi- 


sinusoidal variatiorQ). In addition, iBravo et al.l (|2014ll 
performed an unprecedented comparison of different vari¬ 
ability signatures in wavelet maps, including rotational 
modulation. This approach was also considered in our 
visual inspection procedure for better identification pho¬ 
tometric variability induced by spots. As another se¬ 
lection step, we considered the fact that in certain LCs 
(approximately 5% of a sample) the rotation period may 
correspond to some multiple of the main periodogram 
peak, instead of the ma in peak itself, as explained in 
IDe Medeiros et ^ (|2013fl . This phenomenon may be ob¬ 
served simultaneously in the Lomb-Scargle and wavelet 
periodograms for the same LC. lDe Medeiros et al.l (|2013[1 
proposed that in such cases, the phase diagram of the 
LC typically exhibit two sub-cycles of different depths, 
which compose a full cycle. For the sake of simplicity, 
we excluded such cases through visual inspection. 

From the 408 KCP stars considered in our parent sam¬ 
ple, only LCs with S/N > 1.0 that exhibited clear semi- 
sinusoidal signatures were selected. This selection re¬ 
sulted in a final sample of 131 stars. Finally, we flagged 
LCs for which the periods derived using the Lomb- 
Scargle and wavelet methods agreed with ^10%. These 
were used as an high confidence subsample for evaluation 
of the results. This subsample of high confidence periods 
comprised 67 of the 131 stars (see Table |T|). 


2.2. Pre-selection of KOIs 

Because the KOIs were used as a comparison sam¬ 
ple, a stricter automatic pre-selection was applied. Only 
LCs with S/N > 2.0 for which the periods derived using 
the Lomb-Scargle and wavelet methods (agreed within 
~10% discrepancy) were selected. This automatically 
pre-selected sample was then visually Inspected In the 
same manner as for the KCP host stars. In this step we 
selected only LCs with transits fainter than the peak- 
to-peak amplitude of the semi-sinusoidal variability. Of 
the 3, 798 KOI stars considered In our parent sample, 409 
objects were selected for comparison with our KCP sub¬ 
sample. 


^ In short, semi-sinusoidal variability is characterized by some 
asymmetry of the maximum and minimum fluxes with respect to 
the average flux over time, somewhat irregular long-term ampli¬ 
tude variations, and semi-regular multi-sinusoidal short-term flux 
variations with a typical amplitude of < 0.5 mag and typical period 
of > 0.3 d. This description was developed based on the dynamic 
behavior of star spots, as observed on the Sun. 
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Figure 2. Typical periodograms and wavelet maps for several stars from our final sample. Each panel is configured as follows: at the 
top is the LC (distinguished by its KIC identifiers), at the bottom left is the Lomb—Scargle periodogram; and at the bottom right are the 
wavelet map and the global spectrum. In the Lomb-Scargle periodograms, the power corresponding to a false alarm probability (FAP) of 
0.01 is not displayed because it is visually very close to the x—axis. 
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Table 1 

A sub—sample of the final KCP sample with the rotation periods computed in this work^. 


KIC 

Prot^ 

(d) 

Amplitude'^ 

(e~ cadence"^) 

(K) 

log 9“^ 

(cm s“^) 

(Mq) 

R-star^ 

(R©) 

S/N 

Flag^ 

757450 

19.19 ± 2.39 

0.00383 ± 0.00296 

5152 ± 93 

4.54 ± 0.30 

0.88 

0.88 

8.0 

1 

2165002 

22.59 ± 2.03 

0.00147 ± 0.00214 

5245 ± 83 

4.57 ± 0.30 

0.8 

0.79 

2.89 

1 

2302548 

12.35 ± 2.14 

0.00694 ± 0.00297 

(...) 

(...) 

0.87 

0.79 

34.81 

1 

2438264 

17.93 ± 1.87 

0.00032 ± 0.00108 

5031 ± 74 

4.57 ± 0.30 

0.74 

0.81 

1.52 

1 

2692377 

21.51 ± 1.24 

0.00049 ± 0.00115 

(...) 

(...) 

0.99 

1.11 

4.96 

2 

3323887 

16.81 ± 0.69 

0.00064 ± 0.00132 

(...) 

(...) 

1.07 

1.02 

4.08 

1 

3541946 

15.21 ± 0.86 

0.00016 ± 0.00370 

5725 ± 71 

4.56 ± 0.30 

0.93 

0.94 

1.07 

2 

3656121 

13.64 ± 0.28 

0.00019 ± 0.00069 

6212 ± 98 

4.38 ± 0.30 

1.09 

1.21 

1.17 

2 

3832474 

16.19 ± 1.77 

0.00388 ± 0.00363 

(...) 

(...) 

0.99 

0.95 

8.94 

1 

3967760 

15.11 ± 1.45 

0.00294 ± 0.00297 

5372 ± 113 

4.48 ± 0.30 

0.95 

0.94 

6.18 

2 


^ Full table available in the online version. A portion is shown here for guidance regarding its form and content. 
^ Uncertainties from Monte Carlos Markov chains estimation, see §[lii 
^ Values from fPinsonneault et all ]2012ll . 

^ Values from KIC. 

® Values from KCP Public Archive. 

^ Levei of period reliability (1 — higher, 2 — lower). 


2.3. Angular momentum 

For all of the stars in our final sample the stellar an¬ 
gular mome ntum was estimate d following the procedure 
proposed bv lAlves et al.l (1201011 . When possible, we also 
estimated the total angular momentum, i.e., the angular 
momentum of the host stars plus the contribution of its 
planets, using the following simple central force model; 

Jstar-pl = M \/G{Mstar + Mpi)a{l - 6^) (1) 

where /r is the reduced mass, G is the gravitational con¬ 
stant, Mstar is the mass of the host star, Mpi is the total 
mass of all planets in the system, a is the semi-major 
axis of the orbit, and e is the eccentricity of the orbit. 
The uncertainty in the angular momentum was calcu¬ 
lated via error propagation, using the uncertainties in 
the stellar mass and radius provided in the NASA Exo¬ 
planet Archival. 

2.4. Error estimation 

In principle, we can simply derive Prot from the out¬ 
put of the Lomb-Scargle periodograms and the wavelet 
method, and then use these results to extract other 
features (i.e., the amplitude and p hase) through least 
squar es fitting. Some authors (e.g., iWalkowicz fc Basril 
suggest estimating the period un certainty from 
the periodogram, wher eas others (e.g., iKovacsl 119811 : 
IHorne fc Baliunaall986D suggest using a combination of 
features from the light curve data, a fit of the rotational 
period, and the periodogram. 

We propose a different approach that is based on 
Bayesian analysis, which provides a natural framework 
for the estimation of both a model and its uncertainty. 
The method proceeds as follows: the prior probability 
distribution function provides the method with knowl¬ 
edge of the parameters and their uncertainties before the 
observational data are incorporated. The likelihood func¬ 
tion provides the information regarding the data itself. 
The posterior probability distribution function is then 
constructed from these two inputs. 

® http://exoplanetarchive.ipac.caltech.edu/ 


We introduced the parameters from the period fit (Tj-ot, 
amplitude, phase) into the prior probability distribution 
function. We also introduced the possibility that the 
calculated Prot may be twice or half the real rotational 
period (see § 2.2.2 in iDe Medeiros et al.ll2013l ). The po¬ 
tential errors in the flux, which were calculated from the 
Kepler Flexible Image Transport System (FITS) files and 
nearest-neighbor flux dis persion in the LCs (see Eq. 1 in 
iDe Medeiros et all 120131) . were introduced into the like¬ 
lihood function. 

Markov chain Monte Carlo sampling was used to im¬ 
plements a modificati on of the Metropolis-Hasti ngs al¬ 
gorithm developed by |Goodm^&|^^ari ((JOlQ)- This 
approach offers the advantage of requiring only sim¬ 
ple hand-tuning and thus, producing samples more 
rapidly than traditional algorithms. To achieve this 
implementation we followed t he procedure proposed by 
iForeman-Mackev et al.l (I2013D . For each LC, 1, 000 itera¬ 
tions were performed to generate a virtual sample, based 
on the information obtained from the Lomb-Scargle 
and wavelet analyses along with our degree of uncer¬ 
tainty. Thus the outcome distribution was rooted in 
our primary calculations. The median of the Bayesian 
Prot distribution differed slightly from our Prot values 
(an average difference of 0.16%), thus its standard de¬ 
viation reflected the period uncertainty of the LC fit (see 
Table [T|) . 

3. RESULTS AND DISCUSSION 

We identified a total of 477 stars that were in 
common with p r evious w orks r e lated to our study 
(iMcQuillan et al.l l2013alR 12014 INielsen et al.l l20lA 
IWalkowicz fc Basril 12013 ; IReinhold et al.l 120131) . For 
those stars, a comparison between our periods and those 
preyiously published is presented in Fig. [31 According to 
this figure, our period calculations are in strong agree¬ 
ment with the literature yalues. This agreement can also 
be obseryed in the upper panel of Fig. [31 in which a Gaus¬ 
sian fit to the histogram of Prot, our work Prot, literature 
yields a full width at half maximum (FWHM) of 0.32 d, a 
yalue that is comparable to the lowest Prot uncertainties 
in our work. From Fig. [31 the existence of harmonics 

at Prot, literature — ^Prot, our work and Prot, literature — 
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Figure 3. Comparison of the rotatio n periods computed in o ur 
work with those from the literature _JMc2uiUai^^alJ_[20^a||H 

1LII2'' . .^ 


20141 : INielsen et al.|[^013l : I Reinhold et aO[2(^^ TWalkowTc^T^Basril 

20131') . The grey area indicates the region that represents a discrep¬ 

ancy of 10% or less, and the dashed lines represents 2 Prot aliases. 
Cumulative frequency distributions (CFDs) are presented for the 
horizontal (panel at the top) and vertical (panel to the right) axes, 
in which the dotted lines mark the median values (both approxi¬ 
mately 18 d). 

Prot, our work/ 2 Can be clcarly observed. It is difficult 
to unequivocally determine Prot in such cases if the S/N 
value is low, but because w e set a sufficiently S/N th resh- 
old (according to § 2.2.2 of iDe Medeiros et al.ll2013ll and 
performed visual inspection of all LCs, we have confi¬ 
dence in our results. A more thorough comparison that 
is focused only on these objects is underway. The lower 
panel of Fig. 0] presents the CFDs of the difference in our 
Prot values and those from the literature, for both the 
KCP host stars and KOIs. It can be inferred that the 
harmonic at Prot, literature ~ 2Prot, our work Can be pre¬ 
dominantly attributed to the KCP host stars, the sam¬ 
ple on which we focused more attention. In addition, 
the results of the Kolmogorov-Smirnov test (distance of 
0.117 and probability of 4.4 x 10“^) and the Anderson- 
Darling test (asymptotic probability of 9. 5 x 10~^) per¬ 
formed using the R sof tware package (IR Core Team I 
120131 IScholz fc Zhull201^ strongly indicate that the two 
distributions are drawn from different parent distribu¬ 
tions. 

Based on Fig. 01 the KCP sample would have shorter 
period values in the present work than in the literature, 
in average, thus being more populated at negative val¬ 
ues of the abscissa than the KOI sample. However, this 
apparent discrepancy is actually a bias produced by a 
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Figure 4. Upper panel: Histogram of the difference in the periods 
from our work and those reported in the literature, from Fig. [3] The 
solid black histogram represents KOIs and the solid blue histogram 
represents KCP host stars. A Gaussian fit with the sum of the 
KOIs and the KCP host stars is depicted by the dashed gray curve 
(mean of 0, FWHM of 0.32, and standard deviation of 0.14) to 
contrast the wings of the black and blue histograms. Lower panel: 
CFDs for the difference in the Prot values determined in our work 
and those presented in the literature for KOIs and KCP host stars. 
The color coding is the same as in the upper panel. Although the 
Kolmogorov-Smirnov and Anderson-Darling tests suggest different 
parent distributions, we verified this is a bias effect (see text). 


pie belonging to the shaded area in Fig. 01 the KOI and 
KCP distributions become highly similar; Kolmogorov- 
Smirnov test yields a distance of 0.0214 and a probabil¬ 
ity of 0.97; Anderson-Darling test yields an asymptotic 
probability of 0.99. 

For the total of 477 publicly available periods that are 
in common with our sample, this work offers a cross- 
validation between our calculated periods and the pre¬ 
viously published values. In addition, we provide stellar 
rotation periods of 63 stars for which values were not 
previously reported in the literature. 


work that are an half of those provided in the literature. 

This is an effect of the distinction between the proce- 
dures for measuring Prot and avoiding their aliases (see 
§ 12.11). In fact, the large data sets from the literature 

behavior of starspot dynamics (Ide Freitas et alJ 120101: 
Bonomo &: Lanzall20I2l lSelloll20I3ll. The main oeriod. 
as observed in the wavelet maps and power spectra, is 
typically the rotation period, as is commonly assumed in 

were typically developed from semi-automatic methods. 

the literature (e.g.. IGarcia et al.ll20I4 IMcOuillan et al. 

A semi-automatic method, which was also performed in 

20I3allbl. I20I4 INielsen et al.l 120131: IWalkowicz & Basri 

the present work to build the KOI sample, is subject to 

20131 iReinhold et al.l 201311, except for approximately 5% 

produce false positives, especially due to aliases. In con¬ 
trast, the KCP sample was obtained with a considerably 
more careful treatment, inspection, and selection of the 

LCs and should have less aliases than the KOI sample. 

To avoid such a bias, we can remove the doubtful cases 
by considering that the shaded area in Fig. 0] represent 
the most confident period measurements for both this 
work and literature. Indeed, if we select only the sam- 

of a samole (De Medeiros et al.ll2013l see also §12.11). The 
second period is an harmonic that may be an effect of the 
superposition of active regions and their dynamics. 

To test our method, we analyzed the hourly total so¬ 
lar irradiance (TSI) and spectral solar irradiance (SSI) 
variation in data from the Variability of Solar Irradiance 
and Gravitv Oscillations (VIRGO. lAndersenl 1199111 in- 
strument on Solar and Heliospheric Observatory (SOHO) 
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satellit^B, as in iLanza et al.l (|2004f) . These authors ver¬ 
ified, based on the wavelet and Lomb-Scargle methods, 
that the rotation period of the Sun can be effectively es¬ 
timated from TSI and SSI data in the solar minimum 
phase, except when the solar photometer (SPM) is cen¬ 
tered at 500 nm of farther from the solar minimum. In 
the case of this exception, the rotational modulation is 
masked by the typical relaxation ti me of the active re - 
gions, which is approximately 5-7 d (ILanza et al.ll2004]l . 
When our method was applied to VIRGO SSI green data 
(centered at 500 nm), the wavelet and Lomb-Scargle 
methods yie l ded r esults that were similar to those of 
ILanza et al.l (|2004D . In this manner, we detected a main 
period of approximately 7.3 d (6.53 d for the wavelet 
method and 8.11 d for the Lomb-Scargle method) for 
a very particular case (ie. at 500 nm and farther from 
the solar minimum, cycle 23). Considering the Kepler 
Response Function (which is centered at 600 nm, with 
5 =1200 II™)) which is compatible with the VIRGO SSI 
green data, we expect to correctly detect the rotational 
modulation of Sun-like stars. 

Figure [5] displays the distribution of rotation period in 
the HR diagram for 85 KCP host stars (upper panel), and 
193 KOI stars (bottom panel), where TL a and log 5 were 
obtained from lPinsonneanlt et al.l (j201^ and the Kepler 
Input Catalog (KIC), respectively. Th eoretical tracks 
for so lar metallicity {Z = 0.014) from lEkstrom et al.l 
(|2012li are included. These tracks span a range of 0.8 to 
120 Mq, beginning with zero-age main sequence (ZAMS) 
stars, with an initial rotational velocity of 0.4 times the 
critical velocity, and ending with the early asymptotic 
giant branch (AGB). The model accounts for the mag- 
netic braking law (see lde Freitas &: De Medeiros! l2013aL 
for a nonextensive approach), and sets abundances that 
are compatible with the solar values. As revealed by 
Fig. El the present stellar sample is composed of stars 
with masses that range from those of low-mass M-type 
stars (namely, the coolest stars in the sample with a min¬ 
imum mass of 0.48 Mq), up to 1.53 Mq. These mass esti¬ 
mations were obtained directly from KCF (3 when avail¬ 
able, and from the NASA Exoplanet Archive otherwise. 
The observed range of rotation periods agrees well with 
the position of the stars in the HR diagram. For example, 
there is an important number of low-mass stars, namely 
those with values of less than approximately 4, 870 
K, for which the period values agree with the distribu¬ 
tion of rotation perio d for M stars prev iously reported 
by other authors fe.g.. llrwin et al.ir2011f) . 

One interesting aspect of the present analysis concerns 
the stars with T^s and logg parameters that are sim¬ 
ilar to those of the Sun. We have identified 5 KCP 
and 12 KOI stars with values within 100 K of 
Teff 0 and logg values within 0.10 dex of loggQ, which, 
within the uncertainties, are consistent with the solar 
values. These stars exhibit rotation period that range 
from 7 to 26 d. Note that the rotation period of the 
Sun ranges from 24.47 days at the equator to 33.5 
days at the poles and is equal to 26.09 days on aver¬ 
age (jLanza et al.ll2003ll . The stars that have Teg and 
log 5 values that are similar to the solar values are the 
KCP stars KIC 3541946, KIC 5351250, KIC 9455556, 

® ftp.pmodwrc.ch/pub/data/irradiance/virgo/ 

II* http://kepler.nasa.gov/Mission/discoveries/ 
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Figure 5. HR diagrams for the KCP (upper panel) and KOI 
(bottom panel) samples, overplotted with evolu tionary tracks that 
account for rotation from lEkstrom et al.l 1120121 1 From left to right, 
the tracks displayed correspond to initial masses of: 1.35, 1.25, 
1.1, 1.0, 0.9 and 0.8 Mq. Their colors represent the theoretical 
Trot values, as indicated by the color scale on the right (in units 
of days), and values outside of this range appear in gray. Filled 
circles represent our data, for which the colors correspond to the 
same color scale used for the tracks. The Sun is represented by its 
usual symbol. 

KIC 11389771, KIC 11709244, and the KOI stars 
KIC 5360920, KIC 6226290, KIC 6392727, KIC 6599919, 
KIC 6869184, KIC 7382313, KIC 7870032, KIC 8676148, 
KIC 9839821, KIC 10811496, KIC 11565544, and 
KIC 12644822 . Of particular interest are KIC 9455556, 
KIC 11565544, and KIC 12644822, which have rotation 
periods that are similar to the Sun’s, namely 25.78 d, 

24.93 d, and 23.51 d, respectively. 

The 523 remaining sources exhibit super- and sub¬ 
solar rotation, with periods that range from 2.05 to 

88.94 d (see Table [T]). Th e period distribution of our 
sample agrees with that of IReinhold et al.l (I2013D . with 
a maximum of approximately 10-20 d. 

When the stellar masses are compared with the cor¬ 
responding Pj-ot values, as in Fig. El the expected decay 
trend is observed: faster rotators have higher masses (en- 
gulfment can be invoked for the rapid rotators with sig¬ 
natures of close planets). However, because of the range 
of masses that we investigated we did not observe the 
change in slope that was detected by iMcQuillan et al.l 
(l2013aD . 

Another notable case is that of the star KI C 12735740, 
whose spectroscopic parameters are given bv lWang et al.l 
(|2013D . This target exhibits values of temperature, sur- 
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face gravity, mass, and radius that are similar to the solar 
values, but it has a lowe r Pmt (19.64 ifc 0.61 d), similar 
to the case of e Eridani (|Metealfe et alJ[20l^ . There is 
a notable spectral similarity between KIC 1 2735740 and 
the S un in the Ca II H region (see Fig. 3 of I Wang et al.l 
l2013f) . which leads those authors to conclude that this 
star is chromospherically inactive. In support of this con¬ 
clusion, the star has a very slow proje cted rotation veloc - 
ity, with V sini (1.43 ± 0.78 km s~^. lWang et al.ll2013ll . 
which is similar to t he solar value (1.6 ± 0.3 km s“^, 
iPavlenko et al.l [2012h . The apparent discordance that 
arises from KIC 12735740 exhibiting a smaller Prot than 
that of the Sun but a similar v sin* value, depends on 
the values of the microturbulence and macroturbulence, 
along with the inclination angle. This angle can con¬ 
tribute significantly to the projected velocity. Therefore 
a strict comparison (which is plausible only in the case 
of homoscedasticity) is difficult to perform. To reinforce 
this claim, it sh ould be noted th at e Eridani has a pe¬ 
riod of 11.20 d (jCroll et 120061) and a v sin* value of 
2.45 km s“^ according to iFischer fc Valeiitil (I2005D . but 
this value decreas es to 1.7 km s~^ in the photometric- 
based measure of iCroll et ^ (|2006lj (see jj 5.2 in the 
cited study). Its temperature of 5,156 K (Idem et al.l 
1200411 is sligh tly sub-solar, nev ertheless, its surface grav¬ 
ity (4.57 dex. iClem et al.ir2004fl is consistent with the so¬ 
lar value. Other relevant parameters of e Eridani include 
its radius and mass, wh ich are similar to the solar v alues 
(0.82 Mq and 0.74 Rq iBaines fc Armstrc)ng|l2(312tl . and 
its younger age of 850 Mvr (iDi Folco et al.ll2004ll .There- 
fore, care must be taken in any direct comparison that is 
based solely on v sin i. 

One of the most notable characteristics of the solar sys¬ 
tem is its large amount of angular momentum, which is 
largely associated with Jupiter, which possesses at least 
two orders of magnitude more orbital angular momen¬ 
tum than the spin angular momentum of the Sun it¬ 
self. In this context, we also analyzed the stellar an¬ 
gular momentum, which was obtained using the relation 
Jstar oc (Mstar/M 0 )“, (this is the stellar contribution 
only, with no planetary contribution) of the stars that 
compose the present sample, to examine how their spins 
compare with that of the Sun. For the a uncertainty 
we performed a Markov chain Monte Carlo analysis, in 
which, Kraft’ s relation (jKraftI IT9671) was employed as 
the fit to the model and fluctuations compatible with 
stellar angular momentum uncertainties were then ap¬ 
plied to the fit. A total of 10, 000 iterations were per¬ 
formed and the uncertainty was defined as in § 12.41 Fig¬ 
ure [3 presents the distribution of the stellar angular mo¬ 
mentum J as a function of stellar mass for our sam¬ 
ple of Kepler stars with planets, for which the J val- 
ues were compute d following the procedure proposed by 
lAlves etahl (1201011 . The Sun is represented by its usual 
symbol. The solid line illustrates the best f it of Kraft’s 
relation Jstar OC {Mstar/MQ)^^ (a=4.9±1.4) (IKrafUll967[) 
with our star sample, where the exponent a of the power 
law was treated as a free parameter. Clearly, the distri¬ 
bution of the stellar angular momentum for stars with 
masses of greater than approximately 0.80 Mq follows, 
at least qualitatively, t he law J oc (M/M ^)'^. as previ¬ 
ously demonstrated by lAlves et | 20 in[ l for stars with 
planets detected using the Doppler method. The mass 



Figure 6. Period—mass diagram for the KCP (red stars) and KOI 
(blue circles) samples. 

range of applicability of this relation is as follows: its 
lower end is given by the mass limit of F- and G-type 
stars (approximately 0.8 M 0 ). Conversely, the upper 
mass lim i t (ap proximately 1.25 M 0 ) was set based on 
iKawaled (I1987D . according to whom such stars possess 
thin convection zones. Above this mass regime the rela¬ 
tion between angular momentum and mass obey s a differ¬ 
ent p ower law (the high-mass relation of § 3 in IKawaled 
mia), because of the emergence of a radiative outer zone; 
consequently, the influence of stellar convective transport 
in stellar envelopes rapidly decreases its influence on the 
relation between angular momentum and mass for high- 
mass stars. 

Of the total working sample of 131 KCP host stars, 
the required parameters for the computation of the total 
angular momentum are present in the literature for only 
38 stars. The bottom panel of Fig. [7] displays the distri¬ 
bution of the total angular momentum of the star-planet 
systems of these stars as a function of stellar mass. The 
Sun symbol in this figure corresponds to the total angular 
momentum of the solar system. Interestingly, although 
most of the planet hosting Kepler stars exhibit a trend 
of excess in stellar angular momentum compared with 
that of the Sun, as observed in the upper panel of Fig. [3 
the angular momentum of the star-planet systems ex¬ 
hibits a somewhat different behavior, with a number of 
KCP host stars exhibiting total angular momenta that 
are comparable to that of the solar system and other 
KCP host stars exhibiting a deficit in total angular mo¬ 
mentum, comp ared with the Sun. 

In addition, lAlves etahl (120101) demonstrated that 
stars that host more massive planets tend to have higher 
angular momenta, whereas stars that host less massive 
planets tend to have lower angular momenta. We tested 
th ese hndings , name ly, the results presented in Fig. 9 
of lAlves et al.l (|2010D . using our sample and, in contrast 
with the results of the previous authors, no clear depen¬ 
dence of the stellar angular momentum on the planetary 
mass was observed. 

There are at least two possible explanations for these 
results: (i) the relatively low-number of stars in the KCP 
sample may hinder a proper statistical treatment of the 
results, and (ii) an important number of KCP sources 
may host undetected planets, thus biasing the results. 
The latter explanation suggests the potential for the de¬ 
tection of new planets in the KCP sample. However, 
another possible explanation f or the discrepancy b etween 
the present finding and that of lAlves et al.l (I 2010 D resides 
in the procedure used for the detection of the planets. 
































































10 


Paz-Chinchon et al. 


The sample selected by lAlves etahl (120101) was based on 
radial velocities, whereas the present sample was chosen 
based on eclipse data, in which the handicap of detecting 
planets closer to the stellar host is not present because 
of the typical observation window of the Kepler satellite. 

4. CONCLUSIONS 

We analyzed photometric variations in the current 
sample of Kepler Planetary Host Stars and deter¬ 
mined Prot for a final sample of 131 sources. These 
periods were calculated independently via the Lomb- 
Scargle and wavelet methods, thereby yielding a reli¬ 
able Prot determination. As remarked in § [21 this unified 
method allowed us to consistently distinguish Prot from 
other manifestations of magnetic activity, even when the 
associated amplitude variations were of the same order 
of magnitude. 

Our Prot values are consis tent with those previous ly 
repor t ed in the liter a ture b y iMcQuillan et 'HT (l2r)13alfp . 
1201^, [Nielsen et"^ (I2013D . iWalkowicz fc Basril (120131) . 
and iReinhold et aLT i 20131) . They are also in strong agree¬ 
ment with the theoretical predictions of lEkstrom et ahl 
(1201^. In particular, the agreement between the tracks 
from lEkstrdm et "ni (|2012t) and our data is clear and pro¬ 
vides good experimental support for stellar rotation the¬ 
ory, despite possible different initial conditions between 
the models and observed stars. 

The present analyses has also revealed an interesting 
group of 5 KCP and 12 KOI stars with Teg and log g 
values similar to those of the Sun and rotation period 
ranging from 7 to 26 days. Of particular interest are the 
stars KIC 9455556, KIC 11565544, and KIC 12644822, 
rotating with periods similar to the Sun’s values, namely 
25.78 d, 24.93 d, and 23.51 d, respectively. 

Finally, the stellar angular momenta of our subsample 
of 131 KCP and 193 KOI stars follow Kraft’s relation, 
thus offering an important generalization of this law for 
a particular sample with known planets and photometric 
Prot measurements. Despite this result, no relationship 
between stellar angular momentum and pla netary mass 
was fo und, in contrast with the results of lAlves et al.1 

(IM^ . 

Open questions remain because robust statistics from a 
larger planetary host sample are required. For example, 
to what extent do low Prot values reflect a bias in de¬ 
tectability? How does planetary systems affect the Prot 
value of the central star? A larger sample will also help 
elucidate the behavior of angular momentum with re¬ 
spect to planetary mass, which is a crucial ingredient in 
the modeling of star-companion systems. Special atten¬ 
tion must be payed to the KOI and KIC stars presenting 
Sun’s rotation rate and Teg and log g solar values, for 
which a solid spectroscopic study could show how close 
of the present day Sun are their evolutionary stages. 
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Figure 7. Upper panel: stellar angular momentum as a function 
of stellar mass for our final sample of 131 KCP host stars (blue 
circles), and a sub-sample of 193 KOI stars (black squares). Mid¬ 
dle panel: total angular momentum of the star-planet systems for 
38 (of the total of 131) KCP stars (green circles), overplotted with 
the stellar angular momentum of the KOI sample (black squares). 
Lower panel: stellar angular momentum versus total angular mo¬ 
mentum of the star—planet systems where the color indicates the 
mass of the host star. A clear dependence of the mass on Jstar is 
observed. In all panels, the Sun is represented by its usual symbol. 
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